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Abstract 

Nuclear effective field theory is applied to the effective range expansion of S-wave nucleon-nucleon 
scattering on a discrete lattice. Lattice regularization is demonstrated to yield the effective range 
expansion in the same way as in the usual continuous open space. The relation between the effective 
range parameters and the potential parameters is presented in the limit of a large lattice. 

PACS numbers: 13.75. Cs, 11. lOGh 
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I. INTRODUCTION 



In the last several years nuclear effective field theory (EFT) has been applied extensively 
to low-energy nucleon-nucleon interactions and to few-nucleon systems At energies below 
1 GeV or so, quantum chromodynamics (QCD) reduces to a hadronic theory containing all 
interactions allowed by the symmetries of the theory. At the very lowest energies, the 
interactions are of contact type among nucleons, with arbitrary number of derivatives. For 
the two-nucleon system, the nuclear EFT has established a concrete systematic foundation 
for the traditional description represented by the effective range expansion (ERE), even 
with large S-wave scattering lengths generated by bound or nearly bound states 3 Q • in 
systems with more than two nucleons few-body contact forces are present — a three-body 
force appears already in leading order — and the EFT provides a well-defined, successful 
extension of the ERE P, (|. At higher energies, pions need to be accounted for explicitly 
in the theory. In this case, the EFT goes beyond the ERE even in the two-nucleon system, 
albeit at the cost of a much more complicated renormalization structure 

Once the leading few-body interactions are determined from few-body systems, the main 
goal of the EFT program is to predict the structure of larger nuclei. Before tackling heavy 
nuclei, one would like to be able to predict the properties of infinite nuclear matter. This 
requires a method of solution whose errors are not larger than EFT truncation errors. A 
few years ago it was suggested that this could be achieved by putting nucleons on a spa- 
tial lattice and using Monte Carlo methods to compute the partition function pj. In this 
first, exploratory investigation we considered two-nucleon contact interactions only, with 
parameters adjusted to nuclear matter properties. Subsequently, works have appeared that 

3 3, H U, Q Il3, Q- As yet, however, a full 



extend this approach in various directions 



application of EFT to nuclear matter has not been carried out. In this work we take the 
first step toward this goal by examining the effective range expansion of nucleon-nucleon 
scattering on a discrete lattice. The extension of this work to the determination of thermal 
properties of neutron matter from parameters from few-body physics is currently underway 

As in any field theory, the parameters that appear in the EFT Lagrangian or Hamiltonian 
are not directly observable, since the separation between them and the high-momentum 
components of loops is arbitrary. This separation is the regularization procedure, such as 
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momentum cutoff and dimensional. Because the separation is arbitrary, relations among 
observables should not depend on the regularization scheme. The relation between EFT 
parameters and observables, on the other hand, does depend on the regularization, and some 
regularization schemes are more convenient to apply than others. The program of predicting 
many-body properties from few-body physics requires that the relation between parameters 
and observables be known within the regularization scheme employed in the solution of the 
many-body problem. Placing nucleons in a lattice is a choice of a regularization scheme. 

In this work, we examine lattice regularization for the effective field theory on a discrete 
three-dimensional cubic lattice of a large size. The special aspect associated with the use 
of a lattice is that the nucleons are interacting in a closed space, different from scattering 
of nucleons in the open space. On this issue, the method of Liischer Il6f is well known in 

n 

lattice QCD, and it has been also studied for the nucleon-nucleon interaction [17(, especially 
on the treatment of the scattering lengths larger than the lattice size. These works examine 
effects of finite-volume lattice space in the limit of vanishing lattice spacing, that is, in the 
continuum limit. Here, we focus on effects of the finite lattice spacing in the limit of large 
lattice volume. 

Our objective is to determine the parameters of the two-body interaction on the lattice 
from known phase shifts. We illustrate the method in the case of sufficiently low energies, 
when the phase shifts can be represented by the ERE parameters. The method is similar 
to the continuum case considered in Ref. [2|, and follows a preliminary attempt invol ving 
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one of us (R.S.) several years ago |l8j ]. An earlier, related work can be found in Ref. 
Our results can be applied to the two-nucleon system at low energies, and the interaction 
parameters thus determined are to be used for the many-body Monte Carlo calculation of 
Ref. |15| . In principle, the whole framework could be used at higher energies, densities and 
temperatures, once pion exchange is included explicitly. 

The cutoff or renormalization scale is kept at finite values, and so is the corresponding 
lattice spacing. As the thermal properties should be examined at the thermal or infinite- 
volume limit, the lattice results needed are at the limit of large lattice space. This is 
usually achieved by performing many-body Monte Carlo calculations with various lattice 
sizes, followed by applying the method of finite-size scaling [20]. Upon application of the 
method, there is no need to consider explicit dependence on the lattice size in our potential 
parameters. Determination of the parameters is greatly simplified at the limit of large space 
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size. In fact, we find that the basic algebra is the same as that in free space, apart from the 
use of the reaction (K) matrix instead of the standard scattering (T) matrix. 

The paper is organized as follows. We briefly discuss the K matrix as a description of 
the two-body interaction in a closed space in Sect. |H] In Sect. 11111 the relation between 
the effective range parameters and the potential parameters is obtained by the use of a 
diagrammatic expansion of the K matrix. An alternative derivation by the direct use of the 
wave function is given in App. |XJ The case of a large, discrete lattice is treated in Sect. IIV1 
An elaboration of the mathematical treatment of the Green's function in this case is given 
in App. El A brief discussion of our results in comparison with Luscher's method and some 
other concluding remarks are presented in Sect. |V| 



II. K (REACTION) MATRIX IN CLOSED SPACE 

Our method is based on essentially the same scattering formalism as the well-known 
is. We find that the use of the reaction, or K, matrix (also termed 



the reactance matrix, or R, matrix) |2ll. l22j|. whose language is more familiar to the nuclear 
physics community, greatly simplifies the formalism. 

We consider two particles of mass M interacting through a potential V(r). The wave 
function for the relative motion, tp p (r), satisfies the Schrodinger equation, 

-(V 2 /M)^ p (r) + V(r)Mr) = E p if> p (r) , (1) 

with E p = p 2 /M = p 2 /M (p = \p\). As Eq. (JTJ) is of second order, ip p {r) can be set to 
describe the physical state of interest as a combination of two independent solutions by 
appropriately choosing boundary conditions. The standard choice of boundary condition is 
that the wave function has, apart from the incident plane wave, an outgoing wave with the 
asymptotic form exp(ipr) /r, or (though less popular) an incoming wave with exp(— ipr)/r, 
either one providing the T matrix, the usual scattering amplitude. Another choice is for the 
wave function to have a standing-wave form, a combination of the two asymptotic forms. 
More explicitly, the wave function in the £-th angular momentum state has the asymptotic 
form: 

Ri{pr) -> ji(pr) - - (S e (p) - 1) hf\pr) 
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for the T matrix, or 

Reipr) -> jt(pr) ~ pK f Xp)n e (pr) 

for the K matrix. Here, K^p) = —(i/p)(Se(p) — l)/(Sg(p) + 1) is the K matrix, and Slip) = 
e 2i5 t j g ^ e g ma ^ r j x expressed in terms of the corresponding phase shift 8^. hf\pr) = 

. Note that the first 



je(pr) + ine(pr) is the spherical Bessel function of the third kind 
term of the spherical Bessel function je(pr) in the above equations forms the incident plane 
wave. From R^ipr), the total wave functions is constructed as 

oo 

Mr) = J2( 2£+ l)i e Reipr)Piicos9) , (2) 
e=o 

where Pe{cos8) is the Legendre polynomial with 9 the angle between p and r. 

Clearly, the choices of the outgoing and incoming boundary conditions are unsuited for 
the description of two particles interacting in a closed space. The choice of the standing-wave 
boundary condition can be made by requiring V(r) and ip p {r) to satisfy periodic conditions, 
such as those that make a cubic box of size L x L x L into a torus, 

V(r + nL) = V(r) , 

tp p {r + nL) = tfj p {r) , (3) 

where n is an integer vector with its components covering a set of all integer values. Equation 
(J3J) restricts the allowed values of {p} to be discrete. In fact, the Green's function satisfying 
Eq. (JSJ) with the standing-wave boundary condition is written as 

Gp(p,r-r) = -j £ * , (4) 



p'(^p) p p 
and obeys 

[-V 2 /M-p 2 /M]G P (p,r-r') = -^ £ ^(r)^(r') = -6(r> - r) + <P P (r)<P P (r') , (5) 

p'(^p) 

where <fi p {r) = exp(ip ■ r). Here, p' = 2-nnjL is the undisturbed (by V(r)) momentum, 
while {p} forms a discrete set of eigenmomenta in the closed space, which are determined 
through a decomposition of the above Green's function as elaborated in Ref. [lfij ]. 

In this work, we examine the large-lattice limit by letting L — > oo. In this limit, the 
periodic condition of Eq. Q becomes ineffective and imposes no special condition, p be- 
coming a continuous spectrum bounded by the inverse of the finite lattice spacing. The p' 



sum becomes an integral, 



1 r d3 p' (a 

fj L ~> *> (2^)3 ' ^ 



where p stands for the principal value of the integral, excluding the contribution from p' = p. 
The range of the integration in Eq. (JHJ) is also restricted by the inverse of the finite lattice 
spacing. 

Our method of using the K matrix is equally applicable to both a large, closed space 
and (open) free space. In fact the formalism and basic algebra are the same. The K matrix 

n 

K(p', p) is defined |22j in terms of ip p (r) with the boundary condition (J2J) as 

K(p',p) = J d^{r)V{r)i> p {r) , (7) 
and it satisfies the integral equation 

K(p',p) = V(p',p) + pj ^V{p',p")G(p";p)K{p",p) . (8) 
Here, V(p',p) and G(p',p) are related to V(r) and Gp(p, r — r') as 

V(p',p) = J dh^{r)V{r)<P P {r) , (9) 
{2nf5\p l! -p')pG{p l ] p) = J St j dh'^i^G^r-r'^ir') 

= (2n) 3 6-\p"-p')p^^, (10) 

respectively. Note p implies p ^ p' in this case. 

Equation (JHJ) is the same integral equation that the standard T matrix T(p',p) satisfies 
for scattering in free space, except for the Green's function satisfying the standing-wave 
boundary condition. Because the two equations are of the same structure, the diagrammatic 
expansions generated from them, as expansions in terms of V(p',p), are also of the same 
structure, apart from the presence of the +ip term appearing in the T matrix. This term 
comes from the p' = p contribution that is included in the T-matrix Green's function 
(usually denoted as G^ + \p',p) for the outgoing boundary condition). The +ip term is vital 
for the T matrix to satisfy the unitarity condition, while the term is not present in the K 
matrix, as the K matrix is Hermitian. 

Successive substitution of Eq. into Eq. (JBJ) yields a diagrammatic expansion of the 
on-shell K matrix. For getting the expansion, however, we must regulate Green's function 
and related (momentum-space) integrals, as discussed in Sect. IIHI 



K(p', p) is expanded in angular momentum states, 

k(p',p) = ~ J2( 2i + w ■ p) K t(p>p) » (11) 

i 

where p' and p are unit momentum vectors. The coefficient (— <Ait/M) is introduced so that 
the on-shell Kf(p,p) is expressed in terms of the £-ih phase shift 5t as 

K t (p,p) = K e (p) = - tanS £ (p) . (12) 
V 

Kp(p) is a real function of p 2 that is known to be analytic around p = 0, so it can be written 
as the effective range expansion with a convergence radius of p 2 ~ 1/(2R) 2 in the case of 



scattering from a potential of the range R 2l| . The S-wave expansion relevant to this work 



is 

K \p) =pcot5 (p) = -- + \r Q p 2 + 0(p 4 ) , (13) 
where ao and r are the s-wave scattering length and the effective range, respectively. 



III. RELATION BETWEEN EFFECTIVE-RANGE AND POTENTIAL PARAM- 
ETERS USING K MATRIX 

In this section, we express the effective-range parameters in terms of the potential param- 
eters using the K matrix without specifying the regularization method. As noted in Sect. |H1 
our method of the K matrix is equally applicable to a large, closed space and to (open) free 
space. The algebra is the same except for the details associated with regularization. With- 
out specifying the regularization method, we can then compare our method to the previous 
works based on diagrammatic expansions of the T matrix, which use different regularization 
methods j^Ellil- In order to solidify the comparison, in App. [X]we also show the derivation 
of the same results using the wave function, instead of the diagrammatic expansion, starting 
from the definition of the K matrix, Eq. 0. The explicit case of the lattice regularization 
(for a large, closed space) will be discussed in Sect. IIVI 

We consider the case where the two particles interact through a short-range potential, 
which is expressed in the form of effective field theory, consisting of a combination of 5 3 (r) 
and powers of the nucleon momentum (square) p 2 , 

V{r) = c (A)5 3 (r) - c 2 (A)[V 2 5 3 (r) + 5 3 (r)V 2 ] + . . . , (14) 
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where the parameters Co, C2, etc. depend on the cutoff scale, A. (In the case of a periodic 
condition (jSJ), 5 3 (r) in V(r) of Eq. (JT3j) is to be replaced by a sum over n of <5 3 (r + nL).) 
In momentum space, 



Here we show explicitly only the leading terms in the potential for the case of interest: 
low-energy phenomena dominated by the S-wave interaction. We do not show explicitly 
higher-order terms such as the P-wave term V • V or relativistic corrections proportional 
to p 4 . A more complete discussion of the various terms can be found in Ref. The 
potential (|14|) is generated by removing from the theory other degrees of freedom, whose 
effects are now subsumed in cq, and higher-order counterterms. For example, in the 
case the particles are nucleons, pion-interaction effects can be effectively included in contact 
interactions for \p\ < m 7r /2, where is the pion mass. 

The potential (jl5j) is singular, in the sense that it requires that the problem be regulated. 
For example, an integral of the Green's function of Eq. (|10|1 becomes 



by the use of a multiplicative regulator F(x 2 ), which satisfies lim^oo F(x 2 ) = and 
\im x ^ F(x 2 ) = l. 

For the sake of comparison with other regularization methods, let us take F(x 2 ) to be 
simply an integrable function of x 2 . We then have 



V{p',p) = c + c 2 (p 2 +/) + ... 



(15) 




(16) 



hip, A) 



2tt 2 



M 




M 



L l {K) + P -R{{p/X) 2 ) 



(17) 



2vr 2 



Here, £i(A) is 




(18) 



with 




(19) 



and R(x 2 ) 




(20) 



S 



We expect Q\ = 0(1), but the exact value depends on the regulator. R(x 2 ) is also a 
regulator-dependent function. For a sharp cutoff regulator F(x 2 ) = 8(1 — x)8(x), we have 

0i = 1 and R(x 2 ) = —ln((l - x)/(l + x)) = 1 + -x 2 + . . . 

2x 3 

Another regulated integral that appears is 

d 3 p' p' 2 F(p' 2 /A 2 ) 



/ 2 (P,A) 



(2vr) 3 A p -A f 



.1/ 



2tp 



L 3 (A)+p 2 / (p,A) 



with 

L 3 (A) = 2vr 2 
For the sharp cutoff regulator, 



d 3 p' 

(2tt) 3 



F(p' 2 /A 2 ) = 6 3 A 3 . 
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We can define analogous integrals J 2n , which satisfy recurrence relations, 

d 3 p' p' 2n F(p' 2 /A 2 ) 



hn(p,A) = p J 



(2 



7T 



-En — -En' 



M 



2tt 



E 2 n+l(A) +p 2 / 2n „ 2 (p,A) 



with 



and thus 



l-2n 



J 2n+1 



M 
"2^2 



(A) = 2tt 



(2tt) 3J 

n 2(n+l) 

5> 2 *L 2(n _ i)+1 (A) + P -^R( P 2 /A 2 ^ 



i=0 



A 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



The diagrammatic expansion of K(p, p) is depicted in Fig. ^ In the case of large S-wave 
scattering length a , c (A)J (p, A) is close to unity and diagrams of all orders in c (A) must 
be included. We denote by K\ the sum of the co(A) contributions, 

Kl = Co (A) + co(A)/o(p, A)co(A) + c (A)/ (p, A)c (A)/ (p, A)c (A) + . . . (27) 

On the other hand, c 2 (A) should be treated perturbatively Q|. We denote the sets of 
diagrams with one insertion of c 2 (A) by A 2 ,A 21 , and Am if they have, respectively, no 
Cq (A) factors, Cq(A) factors either before or after the c 2 (A) insertion, and Cq(A) factors both 



9 




K 2 K21 K121 

(b) 



FIG. 1: Diagrammatic representation of Ki, K 2 , K 2 \, and K\i\. The upper figure (a) depicts K\ 
as a sum of the iteration of Co, Eq. I|27j). with a crossing point as cq. The open bubble is the 
regulated Green's function Iq. The lower figure (b) shows K2,K2i, and K\ 2 \ as in Eq. ()29|) . with 
an open circled vertex as K\ and with a dot vertex as c 2 . The shaded bubble is a regulated Green's 
function weighted with vertex momenta, I 2 . 

before and after the 02(A) insertion. The procedure can easily be extended to higher orders. 
We have 

An 

- jjK (p) = K 1 + K 2 + 2-K 21 + K 121 + 0(cf , p 4 ) . (28) 

After some algebra, we find (suppressing the explicit showing of the p and A dependence 
for a while) 

K 2 = 2p 2 c 2 , 
K 21 = K^c^h+p 2 ^) , 

K 121 = 2K 2 lC2 I I 2 . (29) 

Because 

2 • K 21 + K 121 = 2K X -(I 2 K X + c p 2 I ) , 

Co 
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we obtain for K (p), 



-i 



1 



l + 2(c 2 /co)(p 2 + ^i/ 2 ; 

2 



+ 0(c 2 2 ,p 4 



1 M 
~ + 2r? 



c 

1 



+ 



M 



c 2tt 2 



U + 2^L 3 
c 



-2^ 


-p2 


m l1 


c 


.c 


2^ 



M MM 9 C2 



+ 0(c 2 ,p 4 ) 
p 2 + 0(4p 4 ). 



(30) 



We emphasize that in this derivation, 1 — c io> n °t c , is treated perturbatively. We add 
here a note that Eq. (p^Ojl can be also written as 

c + 2c 2 p 2 



4ir 

-—K (p) = 
M 0W 1 - (c I + 2c 2 / 2 ) 



(31) 



with the understanding that 1 — cqIq (but not cq) and c 2 are treated perturbatively. Up to 
the 0(c 2 ) order, Eq. (j3~Tj) has the same structure as that of K (p) obtained [3] by the power 
divergence subtraction (PDS) scheme (with their definition of c 2 being half of ours). 
The S-wave scattering length and effective range are thus expressed as 



M 1 

4:71 ClO 

M 

-To 



c (A) 
c 2 (A) 



+ 



M 
M 1 



Li(A)4 
R(0) + 



2^£.(A)) + ... 

c (A) 



16tt' v c§(A) 4tt 2 A 
where terms up to next-to-leading order are explicitly shown. Equation 



with Ref. Q]. Note that terms beyond this order must include potential terms of p A and 
higher. 

In this case of a large S-wave scattering length, the dimensionless parameter co(A)6 l 1 AM 
is near its (unstable) fixed point — 27r 2 |24j, 



(32) 



is in agreement 



4lX 

C °< A > = M 



1 

flfi 



20iA 



7T 



+ . . . 



2tt 2 



(33) 



9 t MA ' 

while it flows to the trivial fixed point at the zero value for 7r/(2#iao) 3> A — > 
This observation is consistent with the counting rule that we have followed: Since our non- 
relativistic Hamiltonian is a momentum expansion based on power-counting rules Q], it is 
an expansion about the trivial point. As the leading term of the expansion about p 2 = 0, the 
delta-function potential has to be treated nonperturbatively in order to describe the physics 
near the unstable fixed point, away from the trivial one. 
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IV. CLOSED, LARGE, DISCRETE LATTICE 



We now examine the effective range expansion on a closed, large, discrete lattice. We 
consider a spatial lattice that is simply cubic with lattice spacing a and volume L 3 = (aN) 3 , 
in the limit of iV — > oo. 

The coordinate is discretized in units of a, 



an 



(34) 



where n = n^fj is again a vector with integer components Ui along the directions given 
by the Cartesian unit vectors r», % = 1, 2, 3. Since 



|d 3 r^^a 3 , (35) 



we have also 



1 



5 3 (r) - -^, . (36) 
The range of momenta is limited to the first Brillouin zone, 



--<&<-, (37) 
a a 

for each momentum component i. In the limit N — > oo the momentum is continuum in this 
interval. The wave functions in coordinate and momentum spaces are related by 

«/» «/« Wa d rf rf 

^(an) = / / / ^(p)e p , 

J~n/a J-n/a J-ir/a [ An ) 

which is seen to satisfy the required periodicity even for a finite n, Eq. © with Eq. (J34|) . 

When the standard four-point difference formula is used, the kinetic-energy operator is 
expressed on the cubic spatial lattice as 

3 

tr 

Ha rn/a d Pl dp 2 dp 3 ~ in .„ 1 
j/a J —n/a J — tt/o 

That is, as an operator in the spatial space, p 2 is represented as 

p 2 -> -ip(p) = - 2 V(l - cos(a Pi )) . (39) 

1=1 

12 



3 1 

V 2 ^ p (r) -> -V-r{^((n + r i )a)+V((n-r i )a)-2V'(an)} 
^i a 

J-w/aJ-n/aJ-w/a [ Zn ) a 



P{p)/a 2 becomes p 2 in the continuum limit a — > 0, while for a finite a, we have 

1 . . 12 



(40) 



Pip) = occurs only when p = 0, and thus we have no problem of fermion doubling. Note 
that Eq. (|39|) leads to the formal operator expression for the regulator, 



a 2 p' 2 

TV 2 



u\i r /a r /a r /a dp'[dp'idp'i j fav'Y 2 



7T 



tt/cl J — ir/a J —ir/a (^^*) 



'((t) -J^Ei 1 -<-«>] 



(41) 



showing explicitly that a cubic spatial symmetry is now imposed. This suggests the identi- 
fication 



71 A 

- ~ A . 

a 



(42) 



We then write 
I (p,a) 



M 1 



IP 



7T P7T P7T 



dxdydz 



2a (2tt) 3 J_ w J_ w J_ n 3 — a 2 p 2 /2 — (cos x — cos y — cos z) 



M 
2t? 



a/ it \ \ 71 / 



pa\" 



Here. 



Li 



7T 



1 



7T /*7T /*7T 



dxdydz 



71 



.a/ %7ia J_ n J_ n J_ % ^ — (cos x + cos y + cos z) a 
where 6\ is introduced analogously to Eq. (fT%|) . and 



(43) 



(44) 



7? 



pa\ 2 \ 1 



7T / 



16 



7T /*7T /* 7T 



dxdydz 



[3 — (cos x + cos y + cos z)][3 — a 2 p 2 /2 — (cos x + cos y + cos 2)] 

(45) 



The evaluation of these integrals requires some care. Io(p,a) of Eq. (|4*3^) is related to 
Watson's triple integral 2f(z), discussed in App. |Bj through 



J ( Pl a) = --2f 3- — 



(46) 



Watson's integral If (2;) has branch points at z = ±3, so the limit of small ap/jr is delicate. 
In App. |B]we find 

ll(3-e)=A + Be + 0(e 2 ) , (47) 
with A = 0.505462 ■ • ■ and B = 0.0486566 ■ • • . Therefore, 1 is 



6 1 = nA = 1.58796- •■ , 



(48) 
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in agreement with Ref 
function is 



. |3| (where 9\ = 271-/77 with rj = 3.956 •••). In addition, the R 

(49) 



R [x 2 ) = —B + O (x 2 ) = 0.754330 • • • + O (x 2 ) . 

Higher-order integrals can be obtained as in open space. For example, 

3 — (cos x + cos y + cos z) 



hip, a) 



M 1 rT rlT /,T 

' — 7T <J — 7T <J — 7T 

M r C K \ 2r , 

w L * U +pIo{p > a 



dxdydz 



3 — a 2 p 2 /2 — (cos a; + cosy + cos z) 



Here. 



so 



L, - 



7T 



Ana 3 



7T /*7T /»7T 



dxdydz 



2vr 2 



03 - 



7T 



7T 



(50) 



(51) 



(52) 



For comparison, as noted in Sect. a sharp momentum-cutoff regulator gives Q\ = 

1, R{0) = 1, and 9 3 = 1/3. 

Following the same algebra as described in the previous section, we then find the inverse 
of the S-wave K matrix Kq 1 (p) expressed in terms of the potential parameters co(a) and 
02(a), combined with Eq. (JT3j) . as 



M 



1 M 

+ 



Co (a) 27r 2 



-2 



c 2 (a) 



cola 

— + ^ p 2 + 0(p A ) . 
a 2 



M /7T 

Co (a) 27r 2 3 la 



(53) 



Let us first examine the case C2 = 0. In terms of a and 61, the K matrix is expressed in 
the same form as that for the continuum. The scattering length a® is given as 

1 



4vr 2^x 
+ 



(54) 



ao Mco(a) a 

in the 0((ap/7r) ) order in the power counting. The co(a) in this lowest order, c^\a), is 
then 



4 0) (a) 



4tt / 1 



2tt 2 fa 



(55) 



Here, the approximated expression corresponds to that at the fixed point, and is valid when 



a 
a 



< 20i = 3.17591- 



(56) 
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When C2 is included, with Eqs. (JSJ, (J4l>j) . and (jBTj) we obtain 

M1 1 + ^^ +2 , 3 « 3 £44V 



An ao Co (a) 27r 2 \ a Va/ Co(a) 

M c 2 (a) Ma n/ . 

I^o = 373 - 7Z^(°) + • • • (57) 



167T c Q(a) 47T 3 

Co (a) and 02(a) are determined from a and r by inverting Eq. (jBTfl . Up to 0(02) order, we 



have 



167T 

where c °\a) is given by Eq. and 



co(a) = 4 0) (o){l + ^f8t7(cS 0) (a)) 5 

c 2 (a) = ^(c 0) (a)) 2 , (58) 



4a 

r? = 1 + -5— #(0) « 1 , (59) 



7r 2 r 



the last approximation valid when 

^2 



«4S(0) = 3 - 27098 '"- <60) 



We now examine numerically the validity of the expressions for co(a) and 02(a), Eq. (|58j) . 



in the case the particles are nucleons. Both spin sing' 



et and triplet scattering lengths for 



3: 



a 0s = -23.740 ± 0.020 fm and 



S-wave two-nucleon scattering are known to be large 
aot = +5.419 ± 0.007 fm, respectively (for the neutron-proton system). In contrast, the 
spin singlet and triplet effective ranges have more natural sizes: tq s = 2.77 ± 0.05 fm and 
rot = 1.753 ± 0.008 fm, respectively. These values are comparable with the range l/m n of 
the interaction, which sets the expected limit of validity of the ERE, 

V < A = °' 37 - - 57 flirl = Pmax ■ (61) 

For optimal results, the momentum cutoff A ~ 71 /a should be set greater than p max . This 
requires 

I- <7T, (62) 
^0 

so that the lattice spacing a should be less than about Tr/p max = 8.5 — 5.5 fm, which is a lax 
limit. The inequalities (|62j). (|56|). and (JfciOj) can be realized in nuclear systems with a fairly 
wide range of lattice spacings. 
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Using Eq. (J33J), we have 



(0)/ \ 



1 . 7r 2 3 r O 

1+ w 



(63) 



The second term in the square brackets can be large if the lattice spacing is small. As a 
numerical example, let us take a = 2.0 fm, which corresponds to it /a = 1.57 fm" 1 . We find 
Co(a) ~ 1.8 Cq "* (a), that is, the value of co(a) increases by about 80% by the inclusion of the 
momentum-dependent term in the potential. While this is important when going to next 
order in a calculation, it does not imply a failure of the EFT expansion. The parameters 
of the EFT Lagrangian, such as cq(o), are not directly observable. The convergence of the 
expansion is assured as long as Eq. (J52*j) is satisfied. In fact, the ratio 



c 2 (a) 



c (a J 



Pr. 



7] 



n 2 6 3 r 
40?a 



(64) 



is numerically small: For example, for a = 2.0 fm, it amounts to about 0.06. This ratio 
in fact vanishes in the continuum limit. The perturbative treatment of 02(a) seems to 
be reasonable and the effective range expansion is properly described by the momentum- 
dependent potential, Eq. (|14j). 



V. DISCUSSION AND CONCLUSION 



We have expressed the effective range parameters in terms of the potential parameters 
up to p 2 on a large, discrete lattice, basically in the same way as in free space. Equation 
(J53|) relates two sets of the parameters, 

{cq (a) and 02(a)} and {ao and r } (65) 

in the L — > 00 limit. 

In the case of a finite L, the relation is complicated because the momentum spectrum 
{p} depends on L and also on co(a) and 02(a), 

M -> {P [L; c (a), c 2 (a)]} . (66) 

As well as depending on L, {p} is discrete because only a discrete set of standing waves 
satisfying the periodic boundary condition, Eq. (JSJ), can exist in a finite closed space. The 
exact spectrum of {p} must be determined numerically for the given c (a) and c 2 (a). This 
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is the basic procedure of Liischer's method. It involves an elaborate computation to relate 
ao and r$ (or phase shifts at {p}) directly from the energy (mass) spectrum extracted 
from lattice computations (as usually attempted in the case of lattice QCD). Liischer's well- 
known formula to relate ao and the lowest ener gy f or a given L is, for example, a perturbative 
expansion about p = at the limit of a — » [16], and corresponds to the simplest relation 
coming out of Eqs. (|B3j) and (|p^j) . 

Our objective was to obtain the relation For this, we took the closed space to 

be large, by letting L — > oo in the torus space, {p} is then continuous without the limit 
a — > and, for a finite a, is limited to the first Brillouin zone, Eq. (J37|) . Note a subtle, but 
perhaps basic point in this work, concerning the N — > oo limit that we have taken. For 
a simple cubic lattice, rotational invariance is broken by both the ultraviolet cutoff a and 
the infrared cutoff L. For a large value of N, however, the discrete version of the kinetic 
operator P(p)/a 2 , P N (p)/a 2 , becomes 




As iV — > oo, the rotational symmetry is approached in P^{p)/a 2 — > P(p)/a 2 for p — > 
0. Though we approach the infinitely large closed space by maintaining the spatial cubic 
symmetry (by increasing N for each of the three spatial components), the corresponding 
momentum spectrum near p = effectively approaches that of spherical symmetry. In 
other words, by removing the infrared cutoff, only momenta near the ultraviolet cutoff know 
of the breaking of rotational invariance, see Eq. (|41|). The regularization at N — > oo meets 
the requirement of preservation of the proper symmetry for the momenta of interest in the 
low-energy theory. We emphasize, thus, that we take the L = aN — > oo limit, which differs 
from the iV — > oo limit with finite L. 

We can then expand the K matrix in terms of p 2 , Eq. (|53p. and obtain the relation between 
the two sets of parameters in Eq. through a direct comparison of the expansions power 
by power, without carrying out Liischer's elaborate algebra. Here, for the p 2 expansion 
involving c and c 2 , one must be careful so as to meet the power- counting rules associated 
with the application of effective field theory. As described in this way, our method simply 
amounts to the standard effective range expansion (with some caveats). Our method of the 
K matrix is indeed equally applicable to both of a large, closed space and (open) free space 
with the same algebra, as elaborated in Sect. IIHI and App. |XJ 
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In conclusion, following the appropriate counting rules for the S-wave nucleon-nucleon 
interaction, we obtained Eqs. (|58)l. and (jo^j) for a large, simple cubic lattice, where Eqs. 

(pf9"|). and (fo"2|) hold. In principle the same method can be pushed beyond the effective 
range expansion, through the explicit inclusion of pions. The expressions so obtained tell 
us how low-energy two-nucleon data determine the dependence of EFT parameters on the 
lattice spacing, and can be applied to Monte Carlo calculations of many-nucleon systems in 
large lattices. 
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APPENDIX A: DERIVATION OF EFFECTIVE RANGE EXPANSION IN 
LARGE SPACE, USING OF THE WAVE FUNCTION 

We are going to apply the integral equation, 

^„(r) = <f> p (r) + J d 3 r'G P (p, r - r')V(r')^ p (r') , (Al) 

with Gp(p, r — r') of Eq. (pij with Eq. The method followed here is essentially the same 
as that of Ref. [27J • 

In order to clarify the derivation, let us first consider the potential with c 2 = 0: 

V(r) = c (A)5 3 (r) . (A2) 
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With this potential, Eq. J7J) yields gives 

A>,p) = c o (A)0 p (r = O). 
ip p (r = 0) is determined as follows: Equation (jAl|) gives 



■4> P {r) = e ipr + c (A)M^ p (r = 0) • p 



d 3 p' 



(27r) 3 p 2 — p 



,12 ' 



or, 



(A3) 



(A4) 



(A5) 



^ p (r = 0) = 1 + c (A)/ (p, A)^ p {r = 0) , 
by the use of the regulated Green's function Eq. (JTHJ). We thus find 

( r = 0) = .}. T . r . (A6) 

P 1 - co(A)/o(p, A) 

Equation ()A6|) confirms that p (t" = 0) represents an S wave. From Eqs. (Jll)) and (|A3|) 

with Eqs. f|A"Bj) . (fT7jl . and (fTHjl . we obtain 



M 

# (p) = -— c o (A)0 p (r = 0) 

Mc (A) f Mc (A) 



A 



^i+(f)^(b/A) 2 ) 



4tt [ ' 2tt 2 

The effective range expansion of Eq. (JT3*|l relates co(A) to the scattering length ao, 

1 4tt 2^A 

+ — ^ . 



a Mc (A) 7T 

We now consider the potential of Eq. (JHJ). Substituting it into Eq. (JJJ), we obtain 
#(p,p) = [co(A) +c 2 (A)p 2 ]Vo(p,A) -c 2 (A)^ 2 (p,A) , 



(A7) 



(A8) 



(A9) 



where ^(p, A) and ^(p, A) stand for W(r — * 0) an d V 2, p (r — > 0) with a cutoff A, 
respectively. Following the same procedure as the one for Eq. (jA2|) above, we find that 
■00 (p, A) and ^(p, A) satisfy the coupled linear equations, 



[1 - c / - c 2 / 2 ] + c 2 / o 02 = 1 , 
[c I 2 + c 2 h] 0o + [1 - c 2 / 2 ] 2 = -p 2 



(A10) 



(For simplicity, we suppress the p and A dependence in 0's, c's, and J's in the rest of this 
appendix.) Equation (|A10|) yields 



tfj = [(1 - c 2 / 2 ) + c 2 I p 2 ] /Det 

4>2 = - [coh + c 2 h + (1 - cqIq - c 2 I 2 )p 2 ] /Det 



(All) 
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where 

Det = l- c I - 2c 2 I 2 + c\ll - cpoh ■ (A12) 

K (p) is then 

-(An/M)K (p) = [c + c% + 2(c 2 - c 2 2 I 2 )p 2 + c 2 2 I oP 4 ] / Det , (A13) 

which is exact, obtained from Eq. ([AlOp . 

We now impose power counting rules by treating c 2 perturbatively and by expanding 
about 1 — cqIq. We obtain 



--1/ 



1 






1 


Co 




1 




— + 


Co 





C9 

!-( 

Co 



'c V ' \2rr 2 A" y " J " ' c\ 



(A14) 



We thus recover Eq. (j32|) . in agreement with Ref. 



APPENDIX B: WATSON'S TRIPLE INTEGRAL 

We define a function of a complex variable z, T 3 (z), as 

where <fi = (0i, <p 2 , 03) and 

A(0) = COS 0i + COS 02 + COS 03 . 

When \z\ < 3, the integrand has poles at the values of <fi satisfying A(0) = z. These 
poles generate in T${z) two branch points at z — ±3 and a branch cut between the two 
points. Because of this structure on the z complex plane, the (asymptotic) expansions 
about z = ±3 are complicated. The expansion about z = 3 — is what we would like to 
find, and for obtaining the K-matrix expansion, we need to consider the principal value of 
the integral. Note that if we were to naively expand the integrand, we would find that all 
coefficients of (3 — z) n — except for n = — in 

Uz) = g (3 _ , r . JL^ £ £ £ (B2 ) 
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diverge with the degree of the divergence worsening as n increases. A mathematical com- 
plication here is that each term in the expansion of the principal-valued integral has to be 
evaluated numerically. 



Previously, the triple 



\z\ > 3 was studied l2i 



integral 1^{z) at z = 3 + was analyzed by Watson 



30 



31 



and for 



in connection to random walks on lattices and to lattice 



dynamics in condensed matter. The expansion about z = 3 + was found to be |29j, |3 



Uz) = A- -^{z - 3) 1 / 2 -B(z-3) + 1 



{z - 3f 2 + ... 



(B3) 



V2tt 4\/27r 
where A, B, etc. are the coefficients of the integer powers. We have found that these 
coefficients have been quoted sometimes incorrectly in the literature: The first term Is(3) = 
A has been expressed analytically in terms of Gamma functions jy} 3 but in apparent 



disagreement with the correct numerical value |29l . |3 



A = 0.505462- 



(B4) 



We also find the coefficient of the third term to be 



B = 0.0486566- 



(B5) 



instead of the value 0.014625- • ■ quoted in Ref. [29]. In the following and the rest of this 
paper, we use the values that we believe to be correct. 

Analytic continuation of T%{z), Eq. (|B3)1 . from the region \z\ > 3 to the region \z\ < 3, 
above and below the branch cut, yields 1${z) and I^(z), respectively: 



lf{z) =A± -U(3 - zf/ 2 + B{3 -z) T ^^(3 - zf' 2 + 



The Plemelj formula |33j] then gives 

1 



(2tt) 



IP 



7T P7T P7Y 



d 3 (f) 



-{l+{z) + X-(z)) = A + B(3-z) + 0((3 



(B6) 



(B7) 



near z = 3 for z < 3. 

In the follo wing , we sketch the derivation of Eq. (jB3|) because the literature describing 
the derivation 29j is difficult to locate, and also because our value of the coefficient of the 



21 



third term disagrees with the original one quoted in Ref. [29], as noted above. We first write 



• 129]. 



J 3 (3 + e) = / dte-^lKt) 
Jo 



(B8) 



where 



J (t) = - / dxe {cosx)t (B9) 
1 71 

is the modified Bessel function j23|. The t integral in Is(z) can be divided into two integra- 
tions, [0, T] and [T, oo] for a large numerical value of T: 

l 3 (z)=l*(z)+l b 3 (z). (BIO) 

is expanded about z = 3, 

2J(3 + e)=/ dte- 3t I$(t)-e [ dt te~ 3t J 3 (t) + . . . , (Bll) 

and is numerically computed for each term in the e expansion using the closed form of Io(t), 
Eq. (|B9|) . In 2~|(3 + e) we use instead the asymptotic expansion 



I (t) = -U= I I f - ^ 



2vrt V 8t 128t 2 



or, 



with 



3 3t 



Edi 
Ti ' 



(B12) 



(B13) 



w - 1 w _ 3 A 33 281 



128' 



1024' 



Introducing the Incomplete Gamma function 



which satisfies 



= — 4> m {x) e" 

m m 



dt t m e 



m„—xt 



-i/2 (a; J 



du e~ u 



X|(3 + e) can then be written as 



2|(3 + c) 



-i-3/2 



(2tt) 3 / 2 



1 \ ^ 

(2tt) 3 / 2 T 1 /2 2-j j"» 

|^E (2 ,+ 1)Tl -^06 1/2 



+T l/2 



(B14) 
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The coefficients of e° and e 1 are combined from Xg(3 + e) of Eq. (jBllj) and X|(3 + e) of 
Eq. (jB14|) . and are then numerically computed for various large values of T. By examining 
the numerical results, we obtain the asymptotic expansion of Eq. flBl, with coefficients (jEIjl 
and ()B5|) . Note that the terms of half-integer e powers come only from X|(3 + e). 
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